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Abstract 

The use of orthogonal signaling schemes such as time-, frequency-, or code-division multiplexing 
(T-, F-, CDM) in multi-user systems allows for power-efficient simple receivers. It is shown in this paper 
that by using orthogonal signaling on frequency selective fading channels, the cooperative Nash bar- 
gaining (NB)-based precoding games for multi-user systems, which aim at maximizing the information 
rates of all users, are simplified to the corresponding cooperative resource allocation games. The latter 
provides additional practically desired simplifications to transmitter design and significantly reduces the 
overhead during user cooperation. The complexity of the corresponding precoding/resource allocation 
games, however, depends on the constraints imposed on the users. If only spectral mask constraints 
are present, the corresponding cooperative NB problem can be formulated as a convex optimization 
problem and solved efficiently in a distributed manner using dual decomposition based algorithm. 
However, the NB problem is non-convex if total power constraints are also imposed on the users. 
In this case, the complexity associate with finding the NB solution is unacceptably high. Therefore, the 
multi-user systems are categorized into bandwidth- and power-dominant based on a bottleneck resource, 
and different manners of cooperation are developed for each type of systems for the case of two-users. 
Such classification guarantees that the solution obtained in each case is Pareto-optimal and actually can 
be identical to the optimal solution, while the complexity is significantly reduced. Simulation results 
demonstrate the efficiency of the proposed cooperative precoding/resource allocation strategies and the 
reduced complexity of the proposed algorithms. 
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I. Introduction 

In multi-user systems, all users compete for resources and can cause interference to each 
other. This makes it impossible for any user to gain more profit without harming other users. 
The traditional information-theoretic studies of multi-user systems are mainly focused on finding 
the corresponding rate regions and do not advise how to actually achieve the best rates for all 
users simultaneously (see [1]- [4] and the references therein). It is, however, evident that the 
performance of multi-user systems depends on the balance among the users in the competition 
for resources. Moreover, the points in the achievable rate region are not all stable, or even 
feasible if the selfish nature of the users is taken into account. Indeed, it is reasonable to assume 
that all users will compete for the maximum achievable benefits at all times, which may render 
difficulties to the implementation of any prescribed regulations against the selfishness of users. 
For example, although an outcome corresponding to the case when one user is forced to sacrifice 
its performance for the benefit of other users can be theoretically justified, it is hard to make 
sure in practice that the sacrificed user will not deviate from the regulation which is unfair for 
him. 

Recently, game theory has been recognized as an appropriate tool for studying multi-user 
systems [5]- [20]. It studies the actions of decision makers (players, here wireless users) with 
conflicting objectives, and predicts the users' decisions on future actions (strategies) and the 
outcome of the game. If users compete with each other, the existence of "stable" outcomes, 
corresponding to the so-called equilibria, can be analyzed [8]- [10]. On the other hand, if there 
is a voluntary cooperation among users, the extra benefits for all users can be obtained. The 
corresponding games are called cooperative games and one of the most popular approaches 
developed for cooperative games is the Nash bargaining (NB) approach [21]. 

Although the use of cooperative game theory to recourse allocation in multi-user wireless 
systems is a recent research topic, there are some results available. A two-user power allocation 
game on a flat fading channel (FFC) is investigated in [13]. It is argued that certain points in the 
utility space of the game (i.e., the information-theoretic rate region of the multi-user system) are 
not achievable from a game-theoretic perspective. It is also shown that the NB solution based 
on time division multiplexing (TDM) increases the benefits of all users as compared to the non- 
cooperative Nash equilibrium (NE) solution. The study is extended in [14] to iV-user systems with 
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frequency selective fading channels (FSFCs) for the case when only spectral mask constraints 
(SMCs) are imposed on the users. The NB solution is derived based on joint TDM/frequency 
division multiplexing (TDM/FDM) scheme. Unlike the FFC case, the allocation of frequency 
bins becomes a major problem on the FSFC. A more complex resource allocation game on the 
FSFC with only total power constraints (TPCs) limiting the total transmission power of each user 
is considered in [15]. A water- filling based algorithm is proposed to search for the NB solution 
in a two-user version of the game. The proposed algorithm bargains in many different convex 
subspaces of the original utility space and obtains one NB solution in each subspace. Then, the 
NB solution with the largest outcome is selected. However, the TPCs render the complexity of 
the algorithm high. 

One more application area of cooperative game theory is beamforming in multiple-input single- 
output (MISO) systems [16]- [20]. A two-user game on interference channel is investigated in 
[16], where user strategies are defined as the choices of beamforming vectors. The superiority 
of the cooperative NB solution over the non-cooperative NE solution is demonstrated, and some 
special points such as sum-rate and zero-forcing points are shown to be unstable from a game- 
theoretic viewpoint. Kalai-Smorodinsky-type solutions of cooperative beamforming games are 
further derived in [18]. For the games on two-user MISO systems, it is also shown in [19] and 
[20] that any Pareto-optimal point in the game's utility space can be realized through a certain 
balance between competition and cooperation among the users. 

Game theory has been also used for precoding design. A non-cooperative precoding game is 
analyzed in [9] under SMCS and TPCs, in which a multi-user FSFC is considered and the optimal 
precoding matrices are derived based on the NE. It is shown that the matrix-valued precoding 
games boil down to equivalent vector-valued power allocation games, and the resulted precoding 
matrices adopt a diagonal structured The existence and uniqueness of the NE is guaranteed if the 
communication links are sufficiently far away from each other, and the NE is more efficient when 
the interference power is relatively low as compared to the noise power. Although non-cooperative 
games do not coordinate users and, therefore, allow for low-complexity and distributed solutions, 
they often lead to quite inefficient results for all users due to the lack of coordination. 

'The precoding matrices were mistakenly expressed in [9] as a product of the inverse fast Fourier transform (IFFT) matrices 
and power allocation diagonal matrices, while they should be expressed only as power allocation diagonal matrices. 
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In this paper, we develop cooperative NB-based precoding strategies for the multi-user wireless 
systems using cooperative game theory^ The main contribution of this paper is threefold. First, 
it is shown that cooperative precoding games boil down to cooperative resource allocation games 
under orthogonal signaling set up, that is, the TDM-based cooperation among users for FFCs or 
the TDM/FDM-based cooperation for FSFCs. The precoding matrices adopt a strictly diagonal 
structure in these cases. Second, we show that the process of bargaining among users can be 
physically realized in a distributed and efficient manner with very low information overhead 
if only SMCs are imposed on the users. Third, efficient algorithms for the precoding/resourse 
allocation games are developed for the case when both SMCs and TPCs are imposed on the 
users. Although the bargaining problem appears to be non-convex, efficient algorithms are 
designed based on a proposed classification of the multi-user systems into bandwidth- and 
power-dominant. Then, different manners of cooperation are developed for each type of the 
systems. Such classification guarantees that the solution obtained for each type of the systems is 
Pareto-optimal and actually can be identical to the optimal solution. Moreover, the complexity 
is significantly reduced as compared to the complexity required for solving the original problem 
using exact algorithms. 

The rest of this paper is organized as follows. The signal model is introduced and the coop- 
erative precoding/resource allocation game is formulated in Section HO The precoding/resource 
allocation strategies for cooperative games with SMCs are studied in Section Hill Section [IV] 
deals with the two-user games with both SMCs and TPCs. Section [V] demonstrates our sim- 
ulation results. It is followed by Section [VI] that concludes the paper. All proofs for Sec- 
tions UH Hnl and [IV] are summarized in Appendices A, B, and C, respectively. This paper is 
reproducible research [25] and the software needed to generate the numerical results can be 
obtained from www . ece . ualbert a . ca/ ~vorobyov/ ProgNB . z ip. 

II. System model and Precoding/Resourse Allocation Game Formulation 
A. System model 

Consider an M-user wireless system in which all users transmit on the same wideband FSFC 
with channel length L where L depends on the channel delay spread and the signal symbol 

2 Some preliminary results without proofs have been reported in [22], [23], and submitted [24]. 
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duration [27] Assuming block transmission with block length N for all users, the general 
signal model for user i can be written as 

M 

Yi = GiHjjFjSj + Gi H i s, Gn, (1) 

where S; is the N x 1 information symbol block of user i, F$ is the N x N precoding matrix of 
user i, G« is the N x N decoding matrix of user i, Hjj is the N x N channel matrix between 
users j and i, n t is the N x 1 zero-mean additive white Gaussian noise vector with covariance 
£ , {n i n?} = of I, af is the variance of s i? and {-) H , E{-}, and I stand for the Hermitian transpose, 
expectation operation, and identity matrix, respectively. The information symbols are assumed 
to have unit-energy and be uncorrelated to each other and to noise, i.e., E{SiS?} = I and 
E{Siiif} = 0, where denotes the matrix of zeros. 

In order to decompose a wideband FSFC to flat fading frequency bins, orthogonal frequency 
division multiplexing (OFDM) is adopted. Specifically, assuming that the block length N is 
larger than the channel length L, introducing cyclic prefix (CP), and performing IFFT and fast 
Fourier transform (FFT) at the transmitter and receiver sides, respectively, the signal model can 
be written as [9], [26], [27] 

y< = G^Fs, + Gihi (2) 

M 

where fij = Yl ^j'iFj s j + D n « 1S m e N x 1 interference plus noise vector of user i before 
the decoder, D is the FFT matrix, 4?^ is the N x N diagonalized channel matrix between users 
j and i with its kth element being the sampled frequency response of the A;th frequency bin. 
Both the desired communication channel Ha and the interference channels Hjj (Vj, j ^ i) 
are diagonalized due to the CP insertion and the multiplication by matrices D H and D at the 
transmitter and receiver sides, respectively. 

Considering the general case when all users treat the interference as additive noise, the noise 
covariance for user i before the decoder can be expressed as 

M 

K t = J5{Mf } = a]\ + f&jiFjFf&ft. (3) 

3 FFC can be viewed as FSFC with L = 1. 
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Then, the Wiener filter is the optimal capacity-lossless linear receiver [28], [29]. Thus, the 
decoding matrix Gj can be found as 

G, = Ff H£ (HjjFjFf H^f + R^)- 1 (4) 

and the mutual information (information rate) that user i can achieve is expressed as [30] 

R t (F h F_,) = log (det (I + Ff $f Rl^F,)) (5) 

where F_j is the set of the precoding matrices of all users except user i and det(-) denotes the 
determinant. 

In practice, all users attempt to maximize their information rates under certain power con- 
strains. For the case of FSFC, SMCs are usually considered to limit the powers that the users can 
allocate on different frequency bin. These power limits are denoted as p™ ax (fc) (Vi G VL M ^k G 
Vt N ) where Vt M = {1, • • • , M} is the set of user indexes, tt N == {1, • • • , N} is the set of 
frequency bin indexes. Although SMCs also bound the total power by the value ^2 pf ax (/c) 
for user i, such bound may be loose compared to possibly imposed total power limit P^ max . 
Thus, TPCs may also be needed. The aforementioned SMCs and TPCs can be mathematically 
expressed, correspondingly, as 

E{\[F lSl ) k \ 2 } = [F,Ff ] kk < P r x (k), Vz G Q M , VA; G Q n (6) 
E{\\F lSl \\ 2 } = Tr{F. t Ff } < i^ max , Vz G Q M (7) 

where [-}kk and Tr{-} denotes the A;th diagonal element and the trace of a square matrix. 

B. Cooperative Precoding/Resourse Allocation Game 

Considering the wireless users as players, the choices of precoding matrices as user strategies, 
and the corresponding information rates i?j's as user utilities, the game model of the precoding 
problem can be written as 

r = |^M) {Fj | i G fiju}, {i?j(Fj, F_j) | i G ^m} |- (8) 

In the non-cooperative case, when the game players (wireless users) do not collaborate, the NE 
is a stable strategy combination of the game that satisfies 

i?i(Ff , F N *) > F^), VF^, V* G n M (9) 
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where is the precoding strategy of user i in the NE, F ^ is the combination of precoding 
strategies of all users except user i in the NE, and F- stands for any possible precoding strategy 
for user i. 

In the cooperative scenario, all users are willing to cooperate with each other and agree on a 
common principle in sharing the resources. If the users choose the NB approach as a cooperation 
principle, they aim at maximizing the Nash function (NF) defined in the cooperative utility space 
(rate region) as [31] 

5= J] (Ri(Fi,F_i) — R'i) (10) 

where R[ is the information rate (the utility) that user i can achieve in the predefined disagreement 
point which the users will resort to if the cooperation breaks up. 

In the NB game, the users need to specify also a manner of cooperation according to which the 
bargaining is performed. It is required that a particular manner of cooperation results in a convex 
utility space. In the literature, the users are assumed to cooperate with each other using orthogonal 
signaling schemes such as TDM for FFCs and joint TDM/FDM for FSFCs [13]- [15]. It allows 
no interference among the users. The main technical reason for considering orthogonal signaling 
is that the rate region of a general interference channel is yet unknown. Moreover, the use of 
orthogonal signaling allows for power-efficient simple receives, while it is indeed reasonable 
to assume that the users are equipped with simple matched-filter-based receivers. In addition, 
if the users are allowed to interfere with each other, they need to exchange the interference 
information to achieve a desirable performance. It may significantly increase the overhead in the 
system as well as it also significantly complicates the transceiver design. Therefore, orthogonal 
signaling is indeed a reasonable choice which is also adopted here. It is worth mentioning, 
however, that orthogonal signaling may be inefficient when the interference among the users is 
low [32]. In this case, the resulted rate region may be a small subset of the actual rate region. 
However, it is proved in [32] that the cooperative bargaining problem becomes convex even 
without orthogonal signaling when the interference among users is small (as compared to the 
channel noise), which renders the problem simpler in this case. Moreover, the NE solution has a 
satisfactory performance in the low-interference situation, and cooperation may not be the best 
choice in this case considering the price paid for coordinating the users [9], [10]. Therefore, we 
focus on the case of high-interference in which orthogonal signaling schemes are efficient, i.e., 



November 15, 2009 



DRAFT 



8 



the cooperative solutions based on orthogonal signaling achieves better performance than the 
non-cooperative solutions. 

It has been shown in [9] that the non-cooperative precoding game can be simplified to a power 
allocation game under orthogonal signaling. The following theorem shows that the cooperative 
precoding game can be also simplified to a resource allocation game if orthogonal signaling is 
used, i.e., TDM is used in the case of FFCs or joint TDM/FDM is used in the case of FSFCs. 

Theorem 1: If the cooperation among users is based on orthogonal signaling, the precoding 
matrix of each user in the cooperative precoding game © which maximize the NF rfTOl) under 
the constraints © and optionally (0) adopts a strictly diagonal structure. 

Proof: See Appendix A. 

Theorem 1 can be interpreted as follows. In order to maximize the utilities of all users in 
the cooperative precoding game, the users must adjust their precoding matrices to achieve the 
following two tasks: (i) coordinating the utilization of frequency bins (the public resources); 
(ii) allocating powers (the individual user resources) across the frequency bins. Therefore, the 
cooperative precoding game is more complex compared to the non-cooperative precoding game of 
[9] where the game is solved by performing only individual power allocations among frequency 
bins. 

For further developments, two general assumptions need to be made: (i) The channel infor- 
mation of the desired channel is known at both the transmitter and receiver sides of user i 
only; (ii) The TPCs are tight when they are taken into account, i.e., P™ ax < ^ pj iax (k). 

III. Cooperative Precoding/Resource Allocation Games with SMCs 
The following NB precoding/re source allocation problem with only SMCs is considered 

max TT (i^F-O - i$ subject to: [F<Ff] fcfc < # ax (A;), Vz e Q M , VA; e Q N . (11) 

A. Cooperative strategies for two-user game 

The cooperative NB precoding/recourse allocation game ([TT)) is first considered for two-users 
only, i.e., M = 2 and Qm = ^2- Any stable point in the utility region can be selected as a 
disagreement point. Since the NE point given by [9] 

Ff = Vdiag(pr x ), Vz G Q 2 (12) 
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is stable, it can be selected as a disagreement point. In (fT2l) . p™ ax = [p™ ax (l) ; . . . ,p™ ax (iV)] T is 
the N x 1 vector of power limits on different frequency bins, i.e., the spectral mask vector, and 
(•) T and diag(-) stand for the transpose and the operator that forms a square diagonal matrix 
by writing the elements of a vector in the main diagonal, respectively. It can be seen from (fi"2l) 
that each user exploits maximum allowed power on all frequency bins to maximize its rate. 

Knowing the disagreement point, the manner of cooperation between users can be chosen as 
the joint TDM/FDM for FSFCs (see the arguments in the previous section). The joint TDM/FDM 
prescribes that any frequency bin can be used by only one user at any time instant, but it may be 
shared by different users throughout the operation time. The joint TDM/FDM can be implemented 
with low complexity and the corresponding rate region is guaranteed to be convex. 

The following theorem about the structure of the optimal precoding matrices of the two-user 
TDM/FDM cooperative game CCD) on FSFCs is in order. 

Theorem 2: The NB-based precoding/resource allocation optimal strategies for the two- 
user TDM/FDM-based cooperative precoding/resource allocation game ([77]) on the FSFCs are 
obtained through time sharing of at most two sets of diagonal precoding matrices denoted as 
{F\,F\} and {F\,F\}. The following conditions must be satisfied 

e {0, V^)l. Vi, Wk, VI; F[F l 2 = 0, V/; [Tt{¥] - F 2 }) 2 G Vz (13) 

where i G O2, I £ f^, k G Qn, and tyi = {p™ ax (/c) | k G Qn} is the set of power limits for 
user i. 

Proof: See Appendix B. 

Theorem 2 states that the joint TDM/FDM-based cooperation on N frequency bins can be 
realized by the time sharing of two diagonal precoding matrices under SMCs. It can be also 
seen from the proof of Theorem 2 (see Appendix B) that only one frequency bin needs to 



be shared. Denote this frequency bin as k* and assume that [F}] fe * fc * = ^pf ax (k*). Therefore, 



[F 2 ]fc*A,.* = [F^^fe* = and [Fjifc*** = ^/pf ax (k*). Assuming that user 1 shares frequency bin k* 
for a portion of time (0 < a < 1), the information rate of user i eVL 2 can be written as 

Ri = a log (det (l + M2*«*!i!i > ) ) +(1 _ a ) l og ( det (\ + Ml^jM!^ ) . ( i 4 ) 



07 / / V V (Jf 



Therefore, the two-user cooperative NB precoding/recourse allocation game (fTTI) with the joint 
TDM/FDM cooperation scheme can be converted to the problem of finding k* and a and be 
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simplified as follows. The information rate for user i given in (fl4l) is the summation of user i's 
information rates on all frequency bins, and thus, can be rewritten as 

Ri= on(k)Ri(k), V? G n 2 (15) 

where Ri(k) = log(l + |<frjj(/c)| 2 p™ ax (fc)/<7 2 ) is the information rate that user i obtains on fre- 
quency bin k by using it exclusively for all times, and oci(k) is the time portion during which the 
frequency bin k is allocated to user i. Note that in the NB solution of the game, < a,i(k) < 1 
(i G f^) hold only for k = k*. Then, taking the logarithm of the NF, the NB solution can be 
found by solving the following convex optimization problem 

max log(i*i - i?f ) + log(i? 2 - ) 

subject to: < a* (A;) < 1, Vz G Q 2 , VA; G £In 

ai(k) + a 2 (k) < 1, VJfe G ^tv; ^ > #f E , Vi G fi 2 (16) 
where _Rf E is the rate that user z obtains based on the NE solution. It is worth noting that the 
last constraint in (fT6l) guarantees that both users can achieve higher rates than Rf E (i G ^2) 
through the joint TDM/FDM-based cooperation. Otherwise the users resort to the disagreement 
point and the cooperation breaks up. 

B. Cooperative strategies for M-user game 

Unlike the two-user case, where the NB solution of the cooperative precoding/resource alloca- 
tion game can be formulated as a time sharing between two sets of precoding matrices, it is much 
more complex to coordinate the users' precoding matrices in the M-user game. The structure 
used for the two-user game can not be directly applied here, especially when the number of 
users is large. Therefore, to solve the M-user game, we first partition time into time slots each 
of length T to make it easier for the users to perform time sharing. Moreover, considering the 
case when the number of users or the channel states change over time, the time partitioning 
enables a timely update of the bargaining solution as long as time slots are small enough. In this 
case, the cooperative solution can be obtained through the procedure summarized in Table HI In 
the following we focus on the second step of the procedure in Table H 

As an extension of Theorem 2, the following theorem is in order. 
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TABLE I 



Procedure for finding cooperative solution in an M-user game. 



1 . Initialization: the NE solution for the precoding matrices is obtained and the NE point is used as a 
disagreement point. 

2. Computation: the cooperative NB solution for the precoding matrices is calculated. 

3. Implementation: Implement the NB solution for one time slot. If any changes of the number of users or the 
channel states are detected during this time slot, go back to step 1 in the next slot; otherwise, repeat step 3. 



Theorem 3: Precoding matrices corresponding to the NB solution of the TDM/FDM-based 
M-user cooperative game on the FSFCs have the form 

Fi = r l (t) v /dia g (pf iax ), e n M (17) 

where is a diagonal matrix with its kth diagonal element 

imuJ 1 ' [*«■«<*>] (18) 

[ 0, if ti [bi(k),ei{k)] 

with bi(k) and Ci{k) representing, respectively, the starting and ending time moments between 
which frequency bin k is allocated to user i in a time slot [0,T]. The following conditions are 
then satisfied 

Ti{t) = I; TiMTjit) = 0, Vi, j G Q M , j ^ i (19) 

where t G [0, T] is the time instant in a current time slot. 

The proof of Theorem 3 is similar to the proof of Theorem 2 and is omitted here. It is 
worth noting, however, that the difference is that unlike the two-user game in which at most one 
frequency bin needs to be shared, different groups of users may share different frequency bins 
in the M-user game. The first condition in (fT9l) states that no frequency bin should be vacant at 
any time, while the second condition in (fl9l) requests that no frequency bin be used by more than 
one user at any time. Moreover, it is the length of [bi(k), ei(k)], denoted as oti(k) = ei(k) — bi(k), 
rather than the specific values of bi(k) and ei(k), that affects the rates of the users. Once the time 
portions a^k) (Vz G Qm, VA; G Qn) are fixed, the order of using frequency bins is not important 
to the users. Thus, the key problem is to calculate the time portions ai{k) (Vi G VLm, VA; G VLm) 
that user i obtains on a frequency bin k. Mathematically, this problem can be formulated as the 



November 15, 2009 



DRAFT 



12 



following optimization problem 

max V log( J R, - R™ ) 

on(k), i£QM,k£flN 

subject to: < a>i(k) < 1, Vi e ty&f, VA; g £7jv 

^ Qi(fc) < 1, VA; G fiyv; R, > i?f E , Vz G Qm (20) 

where i?j is the sum of information rates that user f obtains on all frequency bins, that is, 

* = E ,„ g a + «) = E Q , w log a + m . (21) 

To avoid a centralized channel estimation and information exchange overhead among users 
on the cooperation stage, a distributed algorithm for solving (|20l) is developed next. 



C. Distributed algorithm for finding the NB solution 

The problem (l20l) is a convex optimization problem with a coupling constraint. Therefore, it 
can be solved in a distributed manner using the dual decomposition method. 
The Lagrange dual problem to (l20l) is given as 



max 



£ logCi* - - £ A(fc) ( £ «,(*) - 1 



subject to: < aj(A;) < 1, Vi G ^m, VA; g 

i2i > i?f E , Vz G fi M ; A(Jfe) > 0, VA; G Sl N (22) 
where A (A;) (VA; G f2/v) are the positive Lagrange multipliers. 

The problem (1221) can be further converted into a two-level optimization problem with the 
following lower level subproblems 



max 



log(i2, - R™) - HkH(k) 



ai(fc), fceHiv 

subject to: < a^k) < 1, VA; G fijv, Ri > R^ (23) 
for each user z G f^M> and the higher level master problem 

min u i( X ) + Yl A ( fc ) sub J ect t0: x ( k ) > 0, VA; G Q N (24) 

iGUjif fcGSZjv 

where E/j(A) is the maximum value of the objective function in ( |23l given A = [A(l),..., 
X(N)] T - 
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The dual problem (I23l)-(l24l) can be solved based on a distributed structure with a coordinator. 
Since the original problem is convex, strong duality holds and the solutions of the dual problem 
((22)) and the original problem (1201) are the same if Slater's condition is satisfied [33]. For our 
specific problem, we have the following result. 

Theorem 4: The Slater's condition is guaranteed to be satisfied for the problem rt20l) as long 
as the NB solution exists. 

Proof: See Appendix B. 

Theorem 4 can be used to further simplify the lower level problem (l23l) . Substituting (f2T|) 
into the objective function of the sub-problem (1231) . the latter can be rewritten as 



max log ( V a,(k)R t (k) - R™ E ) - V X(k) ai (k) 

subject to: < a^k) < 1, Vfc E f2jv; ^ a i( k ) R i( k ) > R ^ ( 25 ) 

where Ri(k) = log(l + |$ii(/c)| 2 pf iax (fc)/cr l 2 ) is the rate on frequency bin k for user i. The lower 
level subproblems are solved distributively by the corresponding users. 

The Hessian of the objective function of the problem (1251) can be written as 

V 2 /i(«i) = " ( E a i( k ) R i( k ) - rr T (26) 

\ken N / 

where r = [^(1), . . . , R^N), R 2 (l), . . . , R 2 (N), . . . , i? M (l), . . . , i? M (iV)] r and a, = [«,(1), 

aj(iV)] T . It is straightforward to see that V 2 /j(cKj) is negative definite since Ri(k) > (Vi G 

VLM^k E fijv). Thus, each Lagrange problem (1251) is guaranteed to be strictly convex and a 

unique solution exists. More importantly, the information required for solving the ith subproblem, 

i.e., Ri(k) and Rf E , is local to user i. 

A coordinator is needed to solve the higher level master problem. Since the overhead of the 

information exchange and the amount of computations for (1241) is insignificant, any user can act 

as a coordinator or all users can serve as coordinators in a round-robin manner. The algorithm 

for solving the dual problem is summarized in Table [III Then, the complexity of finding the 

bargaining solution is determined by the complexity of the lower level subproblems (1251) which 

is 0{N 3 ). 

Note that the coefficients A(fc) (k E CIn) have specific physical meaning. Indeed, X(k) 
represents the risk that cooperation among users breaks up due to a conflict on sharing frequency 
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TABLE II 

Dual decomposition algorithm for NB . 

1. The coordinator initializes A = A° = [A°(l), A°(2), A°(7V)] T and broadcasts it to all users. 

2. Each user solves d25l ) according to the present value of A and transmits its solutions for 
cti(k),k <G £In to the coordinator. 

3. The coordinator updates A according to the gradient of the master problem d24b as A(fc) ~ [A(fc) 
—8(1 — a i(k))]+ (V& G ^n) where (-) + denotes the projection onto non-negative sub-space 

and S is the step length of the algorithm. 

4. If |A(fc) — A(fc)| < £ (Vfc € iljv), stop; otherwise the coordinator broadcasts A and go to step 2. 
Here, £ is the stopping threshold of the algorithm. 

bin k. Thus, in the lower level subproblems, the objective for each user consists of two parts. 
On one hand, a larger cti(k) is preferred to increase the total information rate of user z. On the 
other hand, if cti(k) becomes too large, the cooperation may break up and the utility of user z 
will return to the inferior competitive solution. 
IV. Cooperative Precoding/resource allocation games with SMCs and TPCs 

The following NB precoding/resource allocation problem with both SMCs and TPCs is con- 
sidered 

max J] F-O-iO 

l&lM 

subject to: [F,Ff } kk < p?**^), Vz G Q M , VA; G Q N - Tr{F,Ff } < P™ 3 *, Vz G Q M . (27) 
Unlike the problem (fTT|) considered in the previous section, the diagonal elements of the preced- 
ing matrices F* (Vz G Qm) in (1271) do not necessarily satisfy \Fi\k,k — ^/Pmax(k) when frequency 
bin k is allocated to user i because of the total power constraint. However, if the joint TDM/FDM 
cooperation scheme is used, Theorem 1 applies, and Fj (Vz G Q,m) can be written as 

F< = diag(VpI), Vz G n M . (28) 
Using the same train of arguments as in the previous section, (|27T) can be simplified as 
max V log( J R 4 - 

l&lM 

subject to: 0<oii(k)<l, Vz' G VA; G fijv; ^ «i(A;)<l, VA; G f^; R i >R' i , Vz' G fi A/ 
X] oci(k)pi(k) < Pr\ Vz G ft(fc) < ft max (^ Vz G fi M , VA; G ^ 

(29) 
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where = \pi(l), Pt(2), . . . , Pt(N)] is the power allocation vector for user i, i?- is the dis- 
agreement point for user i, and E4 — ^2 log(l + \&u(k)Fu(k)\ 2 / of) = Yl ck^ (A;)log(l + 
\&ii(k)\ 2 pi(k) / af) is the sum information rate that user i can obtain. 

Unlike the problem (|2Q|) in the previous section, it can be seen that Pj (i G VLm) are also 
optimization variables in (|29l) . Moreover, (|29l ) is non-convex. Indeed, the Hessian matrix of 
fi(ai,Pi) = Yl a i(k)Pi(k) can be written as 



Thus, Hf- (Vi G Qm) are orthogonal matrices, i.e., H^HT. = I (Vi G fi A/ ). The eigenvalues 
of the orthogonal matrices can only be 1 or —1. Moreover, it is known that the summation of 
all eigenvalues of Hf. equals Tr{Hj.} which is zero for any i. Therefore, Hf. (Vz G Qm) must 
have equal number of eigenvalues 1 and —1. The latter means that Hf. (Vi G &m) are indefinite. 
Thus, the constraints a i(k)pi(k) < P™ ax (Vi G Qm) are non-convex and the non-convexity 



of the optimization problem ( |29l ) follows. 

In the following studies, the two-user case is considered and the disagreement point is chosen 
at the origin of the rate region, i.e., R[ = (Vi G fi A ./), instead of the NE point since finding 
the NE solution, in this case, is itself a complicated problem. 

A. Bandwidth-dominant and power-dominant systems 

Finding the TDM/FDM-based NB solution of the problem d29l ) requires joint power and 
frequency bin allocation for each user, and the resulting complexity of the two-user game can be 
unacceptably high. Moreover, the TDM/FDM-based cooperation can be inefficient in some cases 
when TPCs are present. To overcome these problems, we categorize systems into two types and 
deal with each type separately. Toward this end, two definitions need to be given first. 

Definition 1: A point x is Pareto-optimal in space S if and only if y = x for all y satisfying 
y y x in S. 

A Pareto-optimal point corresponds to an efficient allocation of system resources. The NB 
solution is one of the Pareto-optimal points in a utility space (rate region). 

For the two-user cooperative game, there is a well known algorithm for obtaining the 
TDM/FDM-based NB solution if only SMCs are imposed on the users [13], [34]. According 




(30) 
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to this algorithm, the frequency bins are first arranged such that Ri{k)/R 2 {k) > Ri(j)/R 2 (j) 
(Vj, k G £In) if k < j, where Ri(k) is the rate that user i can achieve on frequency bin k by 
using pf ax (k) and allowing no interference from other users. Given any integer k E Qjy, let 

ai(fc) = 1, a 2 (k) = 0, Pl (k) =pT x (k), p 2 (k) =0,Wk< k 

ai(fc) = 0, a 2 (k) = 1, Pl (k) = 0, pa(A;) = p^k), Wk > k 

a y {k) = P, a 2 (k) = 1-/3, Pl (k) = p'Tik), p 2 (k) = pT\k). (31) 
Then the point R = [Ri,R 2 ] is guaranteed to be Pareto-optimal in the rate region for any 
< (3 < 1. Varying k and j3, all Pareto-optimal points can be obtained including the NB 
solution of the game with only SMCs. 

Definition 2: All Pareto-optimal points in a convex space S form the Pareto-boundary of S. 

The NB solution for the two-user cooperative precoding/resource allocation game with only 
SMCs can be found by searching on the Pareto-boundary instead of the entire utility space of 
the game. The algorithm of [13] is based on the principle that frequency bins which are "better" 
for a certain user should be allocated to this user prior to the other frequency bins which are 
"inferior". However, this principle may fail and lead to highly inefficient solutions if TPCs are 
also imposed. 

Consider the following simple example. Assume that there are four frequency bins and [Rx(k), 
R 2 (k)} are [0.5,0.1] for k = 1, [2,1] for k = 2, [1,3] for k = 3, and [0.3,1] for k = 4. 
Also assume that p™ ax = [1, 1, 1, 1] and P™ 3 * = 1.5 for both users. Then according to the 
aforementioned principle, the following resource allocation can be obtained «i(l) = a 2 (A) = 
l,ai(2) = a 2 (3) = 0.5,«i(3) = a x (4) = a 2 (l) = a 2 (2) = and Pl {l) = Pl (2) = p 2 {3) = 
p 2 (4) = l,pi(3) = pi(4) = p 2 {l) = p 2 {2) = 0. Note that the TPCs £ a^p^k) < i^ max (i e 
fi 2 ) are used to derive the TDM/FDM coefficients o; 1 (2) = (P™* - p™ ax (l) ai (l)) /p™ ax (2) = 0.5 
and a 2 (3) = (P 2 max -p 2 nax (4)a 2 (4))/p 2 nax (3) = 0.5. The resulting rates are R x = 1.5 and R 2 = 2.5, 
and the point (1.5, 2.5) is obviously not Pareto-optimal. For example, the strategies according to 
which frequency bin 2 is allocated to user 1 and frequency bin 3 is allocated to user 2 for the 
whole time provide higher rates than the allocation performed according to the aforementioned 
principle. It is because the principle in [13] considers only comparative advantages between the 
users, but not the absolute advantages. 

It follows from the above discussion that the presence of TPCs renders a different bargaining 
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problem since there is a need to coordinate between the power allocation and the frequency bin 
allocation. Therefore, a different approach has to be developed. Toward this end, note that we 
can first consider the solutions for the bargaining game with only SMCs (denoted as game Ql), 
and then add TPCs to the game (denoted as game Q2). 

Observation 1: TPCs do not enlarge the utility space of the game. The Pareto-optimal solutions 
for game Ql are also Pareto-optimal for game Q2 if they are achievable. 

Denote the Pareto-boundary of the TDM/FDM utility space of game Ql as V. Then, the 
following proposition is in order. 

Theorem 5: Assume that the frequency bins are ordered such that R\{k) / B,2{k) > R\{j)/ R2H) 
(k,j G Qn) if k < j. A non-empty subset V can be achieved in game Q2 under both SMCs and 
TPCs if and only if there exist 1 < k < N and < a < 1 such that 

k-l N 

^= > a > — = . (32) 

p max (fc) p™ ax (k) 

Proof: See Appendix C. 

According to (|32|) . all multi-user systems can be categorized into bandwidth- and power- 
dominant. If condition (1321) is satisfied, the system is bandwidth-dominant and the rates of 
both users can increase simultaneously only if new frequency bins are added into the system. 
Otherwise, the system is power-dominant and the rates of both users can increase simultaneously 
only when TPCs of the users are relaxed. 

Observation 2: Beginning as a bandwidth-dominant, a multi-user system gradually changes 
towards a power-dominant as the number of available frequency bins increases. 

B. Bandwidth-dominant systems: TDM/FDM based bargaining 

In the bandwidth-dominant systems, the TDM/FDM-based cooperation is efficient in the sense 
that a non-empty subset V can be achieved in game Q2. Denote the Pareto-boundary of the 
TDM/FDM utility space of Q2 as P 2 - Then, for the bandwidth-dominant systems, the bargaining 
can be restricted in the set V = V2 HP only. The resulted NB solution, denoted as S' NB , can be 
sub-optimal as compared to the optimal solution of the non-convex optimization problem (|29|) . It 
is because the power allocation (which is not the dominant factor in this case) is not optimized 
jointly with frequency bins allocation. 
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Denote the optimal NB solution of game ^2 as S^ B . Also denote the TDM/FDM utility 
spaces of games Q\ and Q2 as IA\ and U2, respectively. Then, the following theorem regarding 
the optimality of S' NB is in order. 

Theorem 6: S' NB = S% if Sf B e V. If S' NB ^ Sf B , then Sf B £ V but S' NB e V, which 
means that S°^ B is not Par eto -optimal in IA\ but S' NB is Par eto -optimal in IA\. 

Proof: See Appendix C. 

Theorem 6 leads to the following two conclusions about the optimality of S' NB in the bandwidth- 
dominant systems: (i) S' NB can be identical to the optimal TDM/FDM based NB solution; 
(ii) S' NB is guaranteed to be Pareto-optimal in IA\ (which is larger than U2) even if the optimal 
NB solution is not Pareto-optimal. 

C. Power-dominant systems: FDM/sampled time sharing-based bargaining 

Let us now consider the case of power-dominant systems. The example given in Subsec- 
tion IIV-AI is, in fact, an example of a power-dominant system. From this example, we can make 
the following observation. 

Observation 3: The use of the maximum allowed power on all allocated frequency bins 
generally results in a non-optimal solution for game Q2. 

To verify this observation, let us denote the set of all frequency bins as B, the subset of 
frequency bins which user 1 occupies using the maximum allowed power as £>" iax , the set of 
frequency bins which user 2 occupies using the maximum allowed power as £>™ ax . Then, user 1 
may improve its rate by water-filling on B — £>™ ax instead of using the maximum allowed power 
on £>f ax , while the rate of user 2 can be kept the same. Here B — £>™ ax denotes the difference 
between sets B and £>™ ax , and the general term water-filling is used to represent the specific 
meaning of finding the solution of the following convex problem 

max Yl + £ iU)Pi(j)) subject to: V Pi (j) = /^ max ; Pi (j) < pf^(j), Vj G tl N 

(33) 

which is a single-user multi-channel power allocation problem with constant Si{j) = \ [^u]jj\ 2 /af 
being a measure of the channel j for user i, which depends on the channel gain and channel 
noise. 
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TABLE III 

Algorithm for the power-dominant systems. 

1. Both users i £ O2 perform the water-filling d331 > on B and obtain two sets of frequency bins Bi (i £ fi2). 
Then, B c = Bi f] B2 is the set of frequency bins under competition. 

2. In the first round of this step, user 1 is allocated the set of frequency bins B\ and user 2 performs water-filling 
on B — Bi. In round j (j goes from 2 to L = \B C \ where \B C \ denotes the cardinality of the set B c ), user 1 
selects a subset, denoted as B 3 S , of j — 1 frequency bins with smallest channel gains from the set B c and performs 
water-filling on B~ B{. Then, user 2 performs water-filling on the remaining frequency bins. After the Lth round, 
L points in the utility space are obtained. 

3. Perform L rounds of the aforementioned step 2 for user 2 starting from the state that user 2 is allocated the 
set of frequency bins B2 and user 1 performs water-filling on B — S2. Obtain other L points in the utility space. 

4. Denote the set of 2L points obtained in steps 2 and 3 as T. Find the Pareto-boundary Vt of St where St 
is the minimum convex space containing T. 

5. Bargain on Vt and obtain the solution S'^ B . 



Observation 3 suggests that the power-dominant games have to be played based on a different 
manner of cooperation from the TDM/FDM. A reasonable choice of the manner of cooperation 
is the FDM/time sharing (TS), which considers time sharing between points corresponding to 
different FDM based frequency bin allocation schemes. Then the power allocation, which is 
the dominant problem in this case, is based on the water- filling problem (1331) . However, the 
complexity of finding the FDM/TS based NB solution is high, especially when the number of 
frequency bins is large. To obtain the FDM/TS-based NB solution, the water-filling should first 
be performed for all 2 N possible frequency bin allocations between the users, and the resulted 
2 N points in the utility space should be recorded. Then the TS is used to obtain a minimum 
convex space containing all these points and the NB solution is derived. The complexity of the 
TS is then 0(4^), which is exponential in the number of frequency bins. 

To reduce the complexity, we consider a simplified version of the FDM/TS, which is the 
FDM/sampled time sharing (STS). The proposed FDM/STS scheme finds the optimal FDM/STS 
based NB solution according to the algorithm described in Table Hill 

Let WJ r% {X) denotes the water-filling operator for user i on the set of frequency bins X. 
It returns the maximum rate that user % can obtain by optimizing its power allocation on X. 
Let also the vector of rates R op< corresponding to the FDM/TS-based NB solution S^ B be 
obtained by time sharing of two points (R° ptl , Rf n ) and (R° pt2 , Rf t2 ) in the utility space 
of game Q2, and the time sharing coefficients are A and 1 — A, respectively, that is, R op * = 
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TABLE IV 

The overall algorithm for the two-user NB game with SMCs and TPCs. 

1. Check the condition d32i >: If it is satisfied, go to step 2, otherwise, go to step 3. 

2. System is bandwidth-dominant: Search on the Pareto-boundary V', and return the solution S' N 

3. System is power-dominant: Derive B\, Bi, and B c . Play the 2L rounds and obtain T and Vr- 
Search on Vt, and return the solution S'^ B . 



(\R° P + (1 — X)Rl P ,\R°2 + (1 — A)i?2 P )■ Denote the sets of frequency bins allocated 
to the users in the points {R opt \ R optl ) and {R opt \ R° pt2 ) as {B opt \ B optl ) and {B opt \ B opt2 ), 
correspondingly. Then, the following theorem is in order. 

Theorem 7: The FDM/STS based NB solution S'lf B obtained using the algorithm in Table EH 
can be identical to the FDM/TS based NB solution S^b- V tne y are not identical, the difference 
d between the logarithm of the NF for S°^ B and the logarithm of the NF for S% B is bounded 

/ / WT\B opn ) \ , / WF\B optl ) \\ 
d < mm log o ~ , log , - • (34) 

Proof: See Appendix C. 

The following conclusions can be drawn regarding S'^ B in the power-dominant systems: 
(i) S'x B is the optimal FDM/STS based NB solution. Thus, it is a Pareto-optimal solution in the 
FDM/STS utility space; (ii) S'^ B can be identical to S°^ B ; (iii) The efficiency of depends 
on the ratios u x and u 2 , where u % = WF(B° pt ) /WF{B - Bt). 

D. The two-user algorithm 

The overall algorithm, which combines both the bandwidth-dominant and power-dominant 
cases, for the two-user cooperative NB game is given in Table [IV] 

In the bandwidth-dominant case, the complexity of searching on V is O(N). In the power- 
dominant case, the complexity of the algorithm in Table [Till is determined by the time sharing 
part, which is 0(L 2 ), i.e., the complexity reduction as compared to 0(4^) for the optimal 
FDM/TS based solution (where the time consumed on water-filling is neglected in both cases) 
is dramatically significant, especially for large N. 
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V. Simulation results 

A. Cooperative precoding/resource allocation games with SMCs 

In the first example, we assume that two users share four available frequency bins. The noise 
power cr 2 is 0.01 for both users on all frequency bins. The channel gains of the desired channels 
$11 and $22 are generated as Rayleigh random variables with mean 1. The channel gains of the 
interference channels <fr 12 and $ 2 i are generated as Rayleigh random variables with means 0.7 
and 0.2, respectively. The elements of the spectral mask vector p max are also Rayleigh random 
variables with mean 1. 

In Fig. [Q the NB solution computed according to Theorem 2 is shown together with the 
NE solution. The boundary of the TDM/FDM rate region is also included in the figure. Fig. [2] 
displays the values of the logarithm of the NF under different TDM/FDM frequency bin allocation 
schemes. In this figure, k is the frequency bin being shared and a is the fraction of time that 
user 1 uses the frequency bin k. It can be seen in Fig.[j]that the NB solution lies on the boundary 
of the TDM/FDM rate region and provides significantly larger rates to both users than the NE 
solution. Moreover, the NB solution is fair to both users. It can be also seen in Fig. [2] that the 
largest value of the logarithm of the NF corresponds to the optimal scheme that provides the 
NB solution. 

In the second example, the distributed algorithm for the M-user game developed in Sec- 
tion IIII-CI is tested. It is assumed that four users share six frequency bins. As in the previous 
example, channel gains of the desired and interference channels are generated as Rayleigh random 
variables with means 1 and 0.2, respectively. The elements of the spectral mask vector p max are 
also Rayleigh random variables with mean 1. The step length 5 = 0.2 (if different values are 
not specified) and stopping threshold £ = 10~ 5 are selected. 

The iterations of the NB process are shown in Fig. [3] The four curves on the upper side of the 
figure show the instantaneous information rates that the corresponding users can achieve, and the 
curve at the bottom shows the corresponding values of the logarithm of the NF. The NB and NE 
solutions and the comparison between them in terms of the percentage of improvement provided 
by the NB solution versus the NE solution are shown in Table |V] for one of the runs. It can be 
seen from Fig. [3] and Table IVl that all users obtain supplementary benefit from cooperation. The 
corresponding final allocation of time portions on each frequency bin for each user is shown in 
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TABLE V 

Comparisons between NE and NB 



User 


NE Solution 


NB solution 


Increased by 


1 


1.1296 


2.2707 


101.02% 


2 


1.4014 


2.4906 


77.72% 


3 


1.2952 


2.3992 


85.24% 


4 


1.6957 


2.4175 


42.56% 



Fig. HI It can be seen that frequency bins 1, 2, 3, and 4 are occupied exclusively by users 3, 4, 
1, and 2, respectively, while frequency bins 5 and 6 are shared by users 1 and 4, and users 2 
and 3, respectively. 

Fig. [5] depicts the effect of the step length on the convergence speed of the algorithm. With the 
step lengths 5 G {0.1, 0.2, 0.3}, the corresponding logarithm of NF is shown in each sub-figure. 
It can be seen that the algorithm is time-efficient with a good choice of the step length. 

B, Cooperative precoding/resourse allocation games with SMCs and TPCs 

Fig. [6] shows the system classification according to Theorem 5 versus the total power limits 
and the number of frequency bins for the two-user system. The total power limits of the 
users p™ x an d pmax are e q ua | an( j var y f rom \ to 51. The number of frequency bins N 
increases from 1 to 256. The desired channel gains are randomly generated using Rayleigh 
distribution with mean 1, and the users do not interfere with each other due to the orthogonal 
signaling assumption. The power limits on different frequency bins pf ax (k) (Vz G Cl 2 , VA; G f2/v) 
are uniformly distributed in the interval [1.8, 2.2]. The frequency bins are sorted such that 
Ri(k)/R2(k) > Ri(j)/R 2 (j) (k,j G £In) if k < j. Following the comparative advantage 
based principle introduced in Section IIV-AI the maximum number of frequency bins hi that 
user i can cover is k\ = {maxti|ti G £In, Y^k=i Pi iax (^) < -Pi} f° r user 1 an d ki = 
{max t 2 1*2 G Qn, Ylk=N-t +iP2 iax (^) — ^2} for user 2. The total normalized bandwidth hi 
(with the bandwidth of each frequency bin normalized to 1) that user i can cover is then bx — k±+ 
( p i ~ EfcLi Pr x (k)) / P r x (h+l) for user 1 and b 2 = k 2 + (p 2 - ELn^+i ^ ax (^)) /pT x (N- 
k 2 ) for user 2. Then the variable r = 1 — (61 + 62) /N stands for the system property characteristic 
according to Theorem 5. The system is bandwidth-dominant if — 1 < r < and is power- 
dominant if < r < 1. It can be seen from the figure that the system changes gradually from 



November 15, 2009 



DRAFT 



23 



bandwidth- to power-dominant when new frequency bins are added into the system, while it 
changes gradually from power- to bandwidth-dominant when the total power limits of the users 
are relaxed. 

In our last example, the power-dominant two-user system is considered^! The number of 
frequency bins varies from 4 to 9 (50 runs for each case). The total power limits of the users are 
set as Pi = 2 (i e fi 2 ) for each user, and the power limits on different frequency bins are set to 
1 + x(k) where x(k) is a uniform random variable in the interval [0.2, 0.25]. It guarantees that 
the system is power-dominant. The channel gains on all frequency bins are randomly generated 
for both users using Rayleigh distribution with mean 1. 

Fig. □ shows the FDM/TS- and FDM/STS-based NB solutions S% and S^ B , respectively, 
in 300 simulation runs. It can be seen in the figure that S'^ B is identical to S°^ B for most of 
the cases. Moreover, although the distance between S% B and S°^ B for some cases may appear 
relatively large in the utility space, the difference between the values of the logarithm of their 
NF are small as shown in Fig. [8] Particularly, Fig. [8] depicts the logarithm of the NF for S°^ B and 
S'n B (denoted as NF op< and NF', respectively) versus the number of frequency bins N when the 
total power limits p. max {% e Cj 2 ) are set to 1.5, 2, or 2.5. Every point in the figure is averaged 
over 50 runs. It can be seen from Fig. [8] that the gap between NF op< and NF' is very small, if it 
is not zero. 

VI. Conclusions 

Cooperative NB-based precoding/resource allocation strategies on FSFCs are studied under 
SMCs and optionally TPCs. First, it is shown that the optimal precoding matrices adopt a 
strictly diagonal structure and the NB-based precoding game is equivalent to a corresponding 
resource allocation game if the users are not allowed to interfere with each other, i.e., orthogonal 
signaling is used. The use of orthogonal signaling is practically important since it significantly 
simplifies transceiver design and allows for significant reduction of the system overhead during 
user cooperation. Second, it is shown that the NB solution of the cooperative precoding/resource 
allocation game with only SMCs can be obtained efficiently in a distributed manner (a simple 

4 Note that for the bandwidth-dominant systems, the algorithm for finding the NB solution inherits the algorithm for finding 
the optimal TDM/FDM based NB solution in the precoding games without TPCs which is already studied above. 
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coordinator is required) with inevitable information exchanges among users. The developed two- 
level user cooperation procedure avoids a large system overhead by enabling users to perform 
most of the computations individually using their local information. Third, it is shown that the 
cooperative NB-based precoding/resource allocation game with both SMCs and TPCs is non- 
convex and finding its optimal solution requires joint optimization of the frequency bins (which is 
the public resource) and each user's transmit power (which is the individual resource) allocations. 
The complexity of finding the optimal solution is unacceptably high in this case. Therefore, it is 
proposed to categorize all multi-user systems into bandwidth- and power-dominant depending on 
the bottleneck resource in the system. For different classes of the systems, the algorithms based 
on different manners of cooperation are developed. While the TDM/FDM based cooperation 
is still efficient for the bandwidth dominant systems, the TDM/STS cooperation is used for 
the power-dominant systems. The above classification of the multi-user systems guarantees that 
the solutions obtained by the algorithms corresponding to each category are Pareto-optimal and 
can be even identical to the optimal solutions, while the algorithm complexity is significantly 
reduced. Simulation results demonstrate the effectiveness of the proposed cooperative solutions 
and their superiority to the NE solutions. 

Appendix A: Proof of Theorem 1 in Section m 

The noise covariance in © is equivalent to R_j = of I when the cooperation among users 
is based on orthogonal signaling such as, for example, TDM for FFCs or joint TDM/FDM for 
FSFCsE Thus, Ri(Fi,¥_i) is simplified to 



The Hadamard's inequality (det(A) < Yli a a f° r a Hermitian positive semidefinite matrix A) 
suggests that the determinant in (1351) is maximized when Fj is diagonal. Moreover, the power 
constraints given in © and © are irrelevant to the non-diagonal elements of Fj. Therefore, the 




(35) 



optimal precoding matrices must be diagonal. 



5 It can be, however, any other orthogonal signaling scheme. 
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Appendix B: Proof of Theorems 2 and 4 in Section HIT1 



Proof of Theorem 2 

As follows from Theorem 1, the optimal and F 2 (I G f2 2 ) are diagonal. The three conditions 
in (fl3l) are based on the fact that the joint TDM/FDM is used. First, consider the FDM part. 
Given any division of the frequency bins, both users will use maximum allowed power on all 
frequency bins allocated to them. Thus, the first condition in (fT3l) follows. 

The second condition in (fT3l) is based on the fact that only one user is allowed on any given 
frequency bin at any time. Thus, there must be a user which allocates zero power on any given 
frequency bin at any time. 

The third condition in (fT3l) is based on the fact that the NB solution can be obtained by 

n 

sharing at most a single frequency bin between both users The proof of the latter fact can be 
given by contradiction using the optimality of the NB solution. 

Assume that the NB solution can be obtained only by sharing two or more frequency bins 
between both users and consider the case when two frequency bins m and n are shared. In the 
sharing scheme, user 2 uses a fraction a\ of the time in frequency bin m and a fraction a 2 
of the time in frequency bin n. Let Ri(m) and Ri(n) be the rates that user i can obtain by 
exclusively using frequency bins m and n, respectively. Without loss of generality, we assume 
that R 2 (m) / Ri(m) > R 2 (n) / R 1 (n) . Then either of the following cases must happen: (i) if 
ai_R 2 (m) + a 2 i? 2 (n) > R2(m), there exists 7 G [0, 1) such that ai_R 2 (m) +a2R 2 (n) = i? 2 (m) + 
7-R 2 (n); (ii) if ai.R 2 (m) + a 2 _R 2 (n) < R 2 (m), there exists 7 G [0,1) such that ai_R 2 (m) + 
a 2 R2(n) = 7-R 2 (yi)- 

Case (i) corresponds to the sharing scheme according to which only frequency bin n is shared 
between both users, and user 2 exploits a fraction 7 of time on frequency bin n. According to 
this new sharing scheme, user 2 obtains the same rate on frequency bins m and n as that in the 
original scheme. Then the rate that user 1 can obtain on frequency bins m and n in the new 
scheme is 



(l-7)fli(n) = 1 



> {l-a l )R 1 {m) + {l-a 2 )Ri{n). 




R 1 (n) = (l-a 1 ) 



R 1 (n)R 2 (m) 
R2{n) 



f (l-a 2 )Ri(n) 



(36) 



6 Note that a similar observation has been made in [14], but here we give a different much simpler proof. 
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The last inequality follows from the assumption that R2{m)/R\{m) > R 2 (n) / Ri(n) . It can be 
seen from (l36l) that the rate that user 1 can obtain using the new sharing scheme is equal to or 
larger than that in the original scheme. This contradicts the assumption that the NB solution can 
be achieved only by sharing of two frequency bins between both users. 

A similar result can be derived for Case (ii). Moreover, when more than two frequency bins 
are shared, the above proof can be used iteratively to obtain the same result. Therefore, the 
optimal solution can be obtained by sharing at most a single frequency bin between both users. 
Thus, F, 1 (i e n 2 ) 

can be obtained by adding/deleting a single diagonal element of (i e Sl 2 ) 
and the third condition in (fT3l follows. ■ 

Proof of Theorem 4 

Since the constraints of the problem (1201) are all linear, the Slater's condition reduces to two 
parts with the first part requiring that the feasible domain of / = ^ log(-Rj — -Rf E ) be open and 

i 

the second part requiring that the feasible domain of the whole problem be non-empty. 

It is straightforward to verify that the first part is satisfied. The second part is equivalent to 
the requirement of the existence of the NB solution. This completes the proof. ■ 

Appendix C: Proof of Theorems 5-7 in Section [TV] 
Proof of Theorem 5 

First note that in game Q\ any resource allocation scheme satisfying (l3"TT) results in a Pareto- 
optimal point in the utility space, and vice versa. Thus, the statement of the theorem is equivalent 
to the statement that (l32l) is the sufficient and necessary condition to guarantee that at least one 

set of {k, (3} satisfies the conditions (f3Tb . To prove the sufficiency, let (1321) is satisfied, k = k, and 

fc-i 

(5 = a in (f3Tb . Then the resulting total powers used by the users are P[ = p™ x (k) + ap™ x (k) 

k=i 

N 

for user 1 and P' 2 = ^ p™ ax (fc) + (1 — <5)p™ ax (/c) for user 2. Using d32l) . it is easy to verify that 

k=k+l 

pmax > p/ an( j pmax > p/ Therefore, the sufficiency is proved. The necessity can be proved 
similarly using contradiction. ■ 

Proof of Theorem 6 

The first part of the theorem follows from the independence on irrelevant alternatives property 
of the NB [35]. This property states that bargaining in a convex subset which contains the NB 
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solution of the original set results in the same NB solution. Thus, it is clear that if S' NB ^ S°^ B , 
then S^b & T 3 ' • Since V is the achievable subset of V in game Q2, it is impossible that S°^ B G V 
and S°^ B ^ V simultaneously. Thus, if S°^ B ^ V' , then S°^ B ^ V as well. This completes the 
proof. ■ 

Proof of Theorem 7 

Let R optl > R opt2 . Then R optl < R° pt2 due to the Pareto-optimality. Let also R 1 = (R\,F%) 
and R 2 = (R\,Rl£) be two points generated in steps 2 and 3 of the algorithm summarized in 
Table Hn] such that R\ > R° ptl and R\ > R° 2 pt2 . Denote the sets of frequency bins allocated to 
the users in the points R 1 and R 2 as [B\,B\) and (Bf,B 2 ), respectively. Recalling that R[ = 
(Vi G Q 2 ), i-e-, the disagreement point is the origin, the difference between the logarithm of the 
NF for R op * and the logarithm of the NF for R 1 can be obtained as 

d x =log (Ai?rV(l-A)i?r 2 )+log(Ai?f 1 + (l-A) J Rf 2 )-log( J R})-log(^) 

= log(A^ + (l-A)^)+log(A^ + (l-A)^ 2 



-Kl / V -^2 -'''2 



TjOptl rjOptl \ / r)Opt2 T)Opt2 

< log ( A %- + (1 " A) + log ( A%- + (1 - A)%- 



r\ v ' y ta v ^2 ^2 

< log ^ + (1 _ A)^ = log = log ( ^f^ ) . (37, 

where the inequalities hold because log {\R optl + (1 — \)R opt2 ) +log (\R 2 ptl + (1 — \)R° 2 pt2 ) — 
logGRj) - log(i^) > 0, R optl > RT t2 ,R° 2 ptl < R°2 Pt \ and R\ > R optl , and the last equality is 
obtained by substituting the notations R opt2 = WT 2 {B° 2 pt2 ) and R\ = WT 2 {B - B{). Here 
1 < j < L stands for the index of the round in which R\ is obtained by user 2. 
Furthermore, using the fact that WF 2 {B - B{) > WT 2 {B -6 2 ),©can be simplified as 

B \WF 2 (B-B 2 )J 

It can be derived in a similar way that the difference d 2 between the logarithm of the NF for 
R op< and the logarithm of the NF for R 2 obeys the following inequality 

d 2 < log 1 1 J- . (39) 

Finally, note that neither R 1 nor R 2 have been assumed to be the rates corresponding to the 
optimal FDM/STS-based NB solution. Indeed, R 1 and R 2 are just two of 2L points generated 
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in steps 2 and 3 of the algorithm summarized in Table [TITl respectively. Thus, the rates corre- 
sponding to the actual solution returned by the algorithm are expected to be superior to the rates 
corresponding to the points R 1 and R 2 , or even equal to the rates in the optimal solution R op *|zl 
Therefore, d < min(<i 1 , d 2 ) and d can be equal to zero. This completes the proof. ■ 
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Fig. 2. The logarithm of the NF under different FDM/TDM frequency bin allocation schemes. 
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Fig. 3. Instantaneous information rates and the corresponding logarithm of NF versus number of iterations. 
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Fig. 4. Allocations of time portions on frequency bins {cti(k)}. 
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Fig. 5. The logarithm of the NF versus number of iterations under different step lengths, 5 € {0.3, 0.2, 0.1}. 
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Fig. 6. System classification versus total power limits and number of frequency bins. 
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Fig. 8. The logarithm of the NF for the FDM/TS and FDM/STS based NB solutions. 
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